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Abstract 

The Sommerfeld factor for arbitrary partial wave processes is derived in the non-relativistic 
limit. The s-wave and p-wave numerical results are presented for the case of Yukawa 
interactions. An approximate analytic expression is also found for the Sommerfeld factor 
of Yukawa interactions with arbitrary partial waves, which is exact in the Coulomb limit. 
It is demonstrated that this result is accurate to within 10% for some common scenarios. 
The non s-wave Sommerfeld effect is determined to be significant, and can allow higher 
partial waves to dominate cross sections. 
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1 Introduction 



In systems that contain interacting degrees of freedom, resonances commonly develop which 
cannot be described by a perturbative approach. This is generally known as the Sommerfeld 
effect pQ, and has been widely studied in the fields of chemistry, nuclear and condensed 
matter physics, and also in the dynamics of mechanical and electrical systems. In particle 
physics, effective theories are constructed to deal with bound states, but the main tool for 
evaluating physical quantities such as cross sections and decay rates is perturbation theory. 
It is then common to neglect non-perturbative contributions. This is usually a reasonable 
approximation for relativistic particles, but for slow moving particles the influence of non- 
perturbative scattering interactions before and/or after the main interaction event becomes 
very significant. 

Recently, the importance of the Sommerfeld effect in dark matter annihilation has been 
pointed out in [3j 0J EJ EJ \7\ [8] for properly calculating thermal relic densities, and deter- 
mining signals for indirect dark matter searches. Other examples where the Sommerfeld effect 
is crucial include threshold production of heavy states at colliders, and partial decay rates 
when the products have large phase space suppression. The purpose of this paper is to review 
the origin of the Sommerfeld effect from a particle physics perspective, and demonstrate its 
significance for some common interactions. 

This work was initially motivated in order to determine the Sommerfeld factor for arbitrary 
partial wave (I 7^ 0) processes, which had been overlooked. It is found that, even in cases 
where higher partial waves can safely be neglected at the perturbative level, these channels 
can become dominant when the non-perturbative physics is properly accounted for. For an 
application of the numerical p-wave results in this paper to relic density calculations, see [9]. 
During preparation of this papeJll, a similar derivation of the general Sommerfeld factor has 
been provided in [TO] ITT]. 

Apart from for pure electromagnetic interactions, it is necessary to use numerical simula- 
tions in order to evaluate the Sommerfeld factor. The computational requirements can then 
be an obstacle to including the Sommerfeld effect in relevant calculations. In this paper, an 
approximate analytic expression is also found for the Sommerfeld factor of Yukawa interac- 

x The numerical Yukawa p-wave results were first presented by the author in a talk at the 46th International 
School of Subnuclear Physics, Erice, Sicily, 29 August - 7 September 2008 http://www.ccsem.infn.it/issp2008 . 
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tions for arbitrary partial waves, which is exact in the Coulomb limit. So far as the author is 
aware this is a new result, and is found to be accurate to within 10% for the most common 
applications. 

The structure of the paper is as follows. SectionOoutlines the derivation of the Sommerfeld 
factor for arbitrary partial wave processes. Section[3]then determines the approximate analytic 
expression for the Sommerfeld factor in the presence of Yukawa interactions. In Section U the 
numerical and analytic results are evaluated and compared. Section [5] presents the conclusions. 

2 Derivation of the Sommerfeld factor 

The non-perturbative physics that leads to the Sommerfeld effect can be thought of as the 
limit of perturbative Feynman diagrams with an infinite number of particle exchanges. For 
a two body incoming or outgoing state, it is useful to consider the non-perturbative 4-point 
vertex function, T, in order to quantify the Sommerfeld effect. In this section, the method 
for determining T is first presented. This function is then used to relate the non-perturbative 
cross section of an arbitrary process to the perturbatively calculated prediction. The ratio of 
these results defines the "Sommerfeld factor." 

The non-perturbative 4-point vertex function is a solution of the Bethe-Salpeter equation: 




Figure 1: Bethe-Salpeter equation in diagrammatic form. 

/j4 
— ^ f G{q) G{p x + p 2 - q) T (1) 

where T are "compact" vertices which do not involve any intermediate state composed solely 
of the scattering particles (ie are described by two-particle irreducible diagrams), and G is 
the non-perturbative propagator. The inhomogeneous integral equation for T can be solved 
to arbitrary accuracy by a Liouville-Neumann series or Fourier methods, given that T and 
the propagator are sufficiently well known from perturbative calculations. However, in the 
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non-relativistic limit, such methods are difficult or sometimes impossible to implement. This 
case is most relevant for some common particle physics situations where the Sommerfeld effect 
is of great significance, and so is now concentrated on. 

Fig [T] can equivalently be expressed as the sum of Feynman diagrams with an increasing 
number of T insertions. In the non-relativistic limit, T usually is dominated by single particle 
exchange. For this reason, this class of diagram is often referred to as a "ladder" , (which can 
be better visualised when the scattering particle propagators in Fig [1] are drawn as straight 
lines) . 

The approximation appropriate for the non-relativistic limit is that of instantaneous inter- 
actions [12] . This requires that the dependence of the momenta time components in the matrix 
element is removed. For the intermediate scattering particles that connect subsequent T in- 
sertions, this prescription puts those propagators on-shell. The diagrams with finite ladders 
are then negligible compared to the infinite ladder diagram. The non-perturbative character 
of the vertex function T is now manifest. In the non-relativistic limit, the Bethe-Salpeter 
equation becomes a homogeneous integral equation. To transform this equation into a more 
familiar form, it is convenient to introduce the function x with the following definition: 

x{P\iV2,Pz,Pi) = G{Pi) G{p 2 ) T(p 1 ,p 2 ;p 3 ,p 4 ) (2) 

From now on, the (',p 3 , P4) arguments will implicitly be assumed in x- They do not influence 
the following steps for solving the Bethe-Salpeter equation, which is now expressed as: 

ix(pi,P2) ~ G (Pi) G(p 2 ) J -j—^ ^{p 1 ,p 2 ;q,p 1 + p 2 - q) x(q,Pi + P2 ~ l) (3) 

The leading contribution to T in the non-relativistic limit is usually from single particle 
exchange. As an example, consider a scalar Yukawa theory with the interaction gm^cp*^^, 
where g is a dimensionless coupling constant. In the case of two (j) particles scattering, the 
leading contribution to T is: 

-g 2 mi 

r - 7z br, ( 4 ) 



(Pl -q) 2 -m 



g 2 m 2 



Ipi ~ q| 2 + m2 o - 



(5) 



where u = (Pi)o — <?o an( ^ ^ rie scattering interaction involves single (p exchange. As mentioned 
before, the instantaneous approximation is required for the non-relativistic limit which sets 
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u) = in T. The perturbative vertex function then only depends on the transferred 3- 
momentum, T = U(p 1 — q). This is a general feature in the non-relativistic limit, but it 
should be noted that the x function still retains q dependence. The following parameters 
and function are now introduced: 

P = (Pi-P 2 )/2 P = (p l +p 2 )/2 x(h,k 2 ) = x(h + k 2 ,h-k 2 ) (6) 

In the centre of mass frame, p = (p Q , p), and P = (m^ + E/2 , 0) where E is the total kinetic 
energy of the system. The Bethe-Salpeter equation with the instantaneous approximation 
applied is then given by: 

7^X4 U(p-q)x(P,q-P) (7) 

= G(P + p) G(P - p )j0- i U(p- q') m l') (8) 
It is now useful to also introduce the Bethe-Salpeter wavefunction: 

4 S (q) = / fpxtt*) (9) 

Using this and integrating eq ([5]) over the p variable, the following equation is obtained for 
the scalar Yukawa case after re-arranging a multiplicative factor: 

- E^j 4 s (p) +J-0^V(p- q') ^s(q') = (10) 

where E is the total kinetic energy of the system, U = —Am^ V, and only the leading term 
in E/m^ is kept. The Bethe-Salpeter wavefunction therefore is a solution of the Schrodinger 
equation (here in integral form), with a potential that accounts for the interactions included 
in r. For external states that have non-zero spin, the same steps can be applied and the 
Schrodinger equation is similarly found in the non-relativistic limit. For further discussion of 
the Bethe-Salpeter equation, see [T21, [T3] . 

The result of including the infinite series of scattering interactions has effectively trans- 
formed the incoming states from plane waves to a composite state described by ip BS . The 
corrections to T generally include annihilation channels. This introduces an imaginary term 
into the Schodinger equation which results in a finite lifetime for the composite state. In scat- 
tering processes, bound states cannot be formed without radiating off energy. The composite 
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state found then generally can not easily be described by an effective theory, and its influence 
is usually neglected. 

For an example of including the non-perturbative scattering prior to the main interaction 
event, consider two body annihilation as shown in Fig [21 




Final States 



Figure 2: Diagram with non-perturbative scattering before annihilation 

The non-perturbative matrix element, A4 wlth laddor , is related to the perturbative result, .M w/o 
according to: 

d 4 k 



M with ladder (px ,p 2 ; {p f }) = J -^—^ M Vi/o {p 3 ,p 4 ;{p f })G(p 3 )G{p 4 )T( Pl ,p 2 ;p 3 ,p 4 ) (11) 

where it is assumed that the Sommerfeld effect in the final states is negligible. This neglected 
feature can easily be included by the same technique, but for simplicity it is ignored here. 
The non-perturbative matrix element can be evaluated once the solution of r in eq (pQ) is 
found. In the non-relativistic limit, the instantaneous approximation should be applied to the 
perturbative matrix element, -M w/o . On performing the k integration for this case, it is seen 
that the Bethe-Salpeter wavefunction then weights the perturbative result as follows: 

/d 3 k 
^3 M w/o (k, -k; { P/ }) V BS (k) (12) 

where the centre of mass frame has been chosen. In a partial wave expansion, the leading terms 
in the non-relativistic limit for -M w/o , with a given angular momentum /, behave proportional 
to k l Y lm in the centre of mass frame, where k here is the magnitude of the 3-momentum of 
either particle , and Y lm is a spherical harmonic function. 

For the s-wave case, it can immediately be seen by evaluating the integral that the non- 
perturbative matrix element differs by a factor of ^ BS (r = 0) relative to the perturbatively 
calculated result, where ip BS is the position space representation of the Bethe-Salpeter wave- 
function. The Sommerfeld factor for a given partial wave, I, is defined as: 

non-perturbative partial wave cross section 

Oi = ; ; — ; (t'J J 

perturbative partial wave cross section 
5 



In the non-relativistic limit, S is simply \i^ bs (tc = 0)| 2 . In order to evaluate eq (I12|) for higher 
partial waves, it is convenient to decompose the wavefunction in the (orthogonal) spherical 
harmonic basis. The wavefunction components can then also be written in a separable form 
as defined below: 



^ im (k) = F El (k)Y lm (n k ) 



(14) 



For a particular partial wave process, using A4 w/o = a t k l Y lra in eq (fT2j) with a t being some 
momentum independent factor, the non-perturbative matrix element is equivalently given by: 

dk 



M 



with ladder 



a i J2 

I'm' 



a, 



(0 _ )3 k l+2 F EV (k) J ^ fc n m (^)y; m ,(^)y zw (^ p ) (is) 

/dk 
— 2 k l+2 F El ,(k) 5 lv 5 mm , Y Vm ,{Sl p ) (16) 

r.i„ 



where Q p is the solid angle aligned with respect to the original collision axis of the incoming 
states. For a given partial wave, the Bethe-Salpeter wavefunction generated by scattering 
interactions is forced to have the same quantum numbers. In the perturbative case, the 
momentum appearing in the matrix element is fixed by the incoming states, k = |p 1 | = |p 2 |. 
The Sommerfeld factor can then be determined from eq (I16p using eq (j!3f) to obtain: 

2 



J J* k»* F El {k) | Pl |-< 



(17) 



In order to evaluate the result of integrating the momentum wavefunction weighted by different 
powers of k, it is useful to consider the radial Fourier transformation. This is now determined 
from the usual Fourier transformation: 



^EZm(r) 



d 3 k 

(2^)3 <W k ) e " k ' r 
k 2 dk C 

-— 3 / dQ k F El (k)Y lm (Q k ) J^i-i) 1 3Akr)Y llmt (n r )Yf ml (n k ) 

^ 71 ' J I'm' 

R El (r) Y lm (Q r ) (18) 



The result of differentiation of the radial wavefunction leads to the relation: 



d n R El (r) 

Q r n 



dk 



d n j l (kr) 



(19) 
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A series expansion of the Bessel function then gives, jAkr) = (kr) 1 /(2l + 1)!! + 0(kr) l+2 where 
(21 + 1)!! = (21 + l)!/(2 / 11), and so the following result is found: 



d l R E i( 



Q r l 



r=0 



(2Z + 1)!! 



(20) 



Using this in eq (|17|) . the Sommerfeld factor for the partial wave cross section with angular 
momentum I is equivalently given by: 

2 



(2Z + 1)!! d l R El (r) 



g r l 



r=0 



(21) 



IPil' « 

where | p 1 | is the magnitude of the 3-momentum of either incoming particle in the centre of 
mass frame at infinite separation. This result has been derived in the non-relativistic limit. 
So far as the author is aware, this is a new result that, during preparation of this paper, has 
also been independently presented by Iengo [TU] (and see footnote 1). 



3 Analytic solutions for wavefunctions 

In the previous section, it was shown that the Sommerfeld factor is related to solutions of 
the Schodinger equation in the non-relativistic limit. The relevant potential is constructed by 
considering the interactions present in two-particle irreducible diagrams. For some common 
examples such as gauge boson, meson or Higgs interactions, a Yukawa potential is found at 
leading order. For a non-zero mass of the mediator responsible for the Yukawa interaction, 
there is no analytic solution for the wavefunction. This can then be determined numerically, 
however for many applications, an approximate analytic solution may be adequate for the 
level of precision desired. In this section, an approximate analytic wavefunction is found for 
the Yukawa potential and the Sommerfeld factor is then determined. The results are exact in 
the Coulomb limit, relevant for electromagnetic interactions. 

Given a spherically symmetric potential, V(r), the Schrodinger equation is separable. 
For a two-body system containing particles with a common mass M, the radial part of the 
wavefunction obeys the following: 

(-^ d2 r- M ? 2 + nr) + ^^l)rR l (r) = (22) 

where (3 is the speed of each particle when at infinite separation in the centre of mass frame. 
Natural units (h = 1) are now used for the rest of this discussion. For Yukawa interactions, 
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the potential is of the form: 



Vy = -— (23) 

r 



where is the mass of the mediator for the interaction, and A is the interaction strength 
which is positive (negative) for an attractive (repulsive) interaction. In order to find an 
approximate analytic solution, the following approximation is applied to the potential: 

A 5 e~ Sr 

Vy ~ V H = Y - F ^ (24) 

This approximation maintains the same short and long distance behaviour of the Yukawa 
potential, and is commonly known as the Hulthen potential. The question is then raised as to 
what choice of 8 best reproduces the Yukawa potential. In order to answer this question for 
small strengths of the potential, parametrised by A, it is useful to consider the Schrodinger 
equation in an equivalent form known as the Lippmann-Schwinger equation: 

/•oo 

R t (r) oc j t (kr)+Mk / r' 2 j^kr^ n z (Ar>) V(r') R t (r') dr (25) 

Jo 

fc.r-^O (kr) 1 



M f°° 



(2Z + 1)!! f 21 + 1 J Q ' ' v ' '\ (26) 
where r < = min(r, r'), r > = max(r, r') and k = M(3. The Sommerfeld factor is sensitive to the 
behaviour of the wavefunction at the position space origin. In order for the approximation to 
be precise up to 0{A 1 ) in the limit of zero kinetic energy, it is necessary for the first moment 
of the radial potential, J °° r' V(r') dr' to be unchanged by the substitution V Y — ► V H . This 
fixes the relation: 

7r 2 m, 

6 = (27) 

For k 0, there is no analytic relation between 5 and which equates the first order 
correction to the wavefunction, but limited progress can be made by performing a series 
expansion in k/m+. In the following analysis, the identification in eq (I27p is assumed for 
all strengths of the potential and kinetic energies. However, it should be stressed that more 
accurate results can be obtained if the choice of 5 is improved for these cases. 

The s-wave wavefunctions can now be found analytically, but for I ^ 0, a further approx- 
imation [J2] must be applied to the centrifugal term in order to permit analytic solutions: 

r, - "' - ^ fc.fr «i (28) 

' M (1 — e~ Sr ) 2 Mr 2 y ' 



This is reasonable for short range potentials, but does not reproduce the correct long distance 
behaviour. These approximations convert the effective potential of the radial Schrodinger 
equation into a form that allows solutions in terms of hypergeometric functions, and therefore 
is a special case of the Natanzon hypergeometric potential [15] . For the chosen approximations, 
V H + V, is in the form also known as a Manning-Rosen or Eckart potential. 

The physics of the Schrodinger equation is only sensitive to certain combinations of the 
free parameters. The following dimensionless variables are then constructed: 

w = -— * x = — y = z = 2rMp 29 

x o p m+ 

Introducing a further variable, t = 1 — e~ z ^ 2w , the radial Schodinger equation with the ap- 
proximated Yukawa potential and centrifugal term takes the form: 

The regular solution of this differential equation is given by: 



Ri 



t l+l e -iz/2 



r(A) 



r(o-)r( a +) 



T(2iw) 



2 F 1 {a',a + -\ ] t) (31) 



a ± = l + l + iw(l± y/1 - x/wj A = 21 + 2 (32) 

where 2 F l is a hypergeometric function which has the following behaviour at small t: 

2 F 1 {a+,a-; 2A; t) = l + ^-t + 0{t 2 ) 

The normalisation of the wavefunction has been chosen such that an incident wave e jkr of 
unit strength is obtained at infinity. This is then consistent with cross section calculations. 
In the spherical harmonic basis, an incident wave has the following decomposition: 



e ik.r _r-+oo \ ~« 



oo - I 



elkr _ (.^Ig-ttr] J- Y^r) Ylm(n k ) (33) 



' 2ikr - 

1=0 m=-l 

The limiting behaviour of the 2 F 1 hypergeometric function is given below: 

■p ( — + \ +\ ™ r(A)r(A-a--q+) tz f T(X)T(g- + q+ - A) \ 

Ma > a ;A;t) > r(A-«~)r ( A-a + ) +e { r(a-)r(a+) J (34) 

and so the constructed solution in eq (I3ip has the correct normalisation up to a phase factocl, 
which is irrelevant for determining the Sommerfeld factor. In the Coulomb limit (w — ► oo), 



2 Note the behaviour of the T function under complex conjugation; T(h*) = [r(/i)]* 
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the wavefunction reduces to: 



R l [w — ► oo] 



l p-iz/2 



z e 



r(A) 



?,.r 



(35) 



This result follows from eq (|31|) by considering the limiting behaviour of the T function: 



\T(a + ib)\ 



b^rOO 



'2vr 



» b a e-* b l 2 J— [l + O^- 1 )] fora,6e 



(36) 



where R is the set of real numbers, and the confluent hypergeometric function, is obtained 
using the identity: 



iFxia^b-jv) 



lim 

a n — >oo 



2-Fi(oi,o 2 ; 6; v/a 2 ) 



(37) 



The normalisation of the Coulomb wavefunction is inherited from eq (|31|) , but a check of the 
limiting behaviour of the confluent hypergeometric function confirms that an incident wave 
of unit strength is found as r — > oo. 

The Sommerfeld factor in this dimensionless formalism is given by: 

2 



s, 



1,0 



(2Z + 1)! tfR^z) 



(I 



z=0 



(38) 



However, a naive application of the I ^ approximate Yukawa wavefunction in this formula 
can give results that violate the unitarity limit. This is a consequence of the approximation 
applied to the effective centrifugal potential. A more careful analysis is provided in the next 
section that properly accounts for the centrifugal approximation. 



3.1 Sommerfeld factor with a modified centrifugal term 

After modification of the centrifugal term, the free particle eigenstates are no longer plane wave 
solutions. An equivalent statement is that the canonical momentum of the position co-ordinate 
is modified by the approximation. The transformation between the position representation 
and canonical momentum representation is determined from the usual completeness relation: 

Ir> = /||r Ik')(k'|r) (39) 

Since the spherical harmonics are unaffected by the approximation, the same steps as in 
eq (|18p can be applied to obtain the modified radial transformation relation: 

Rl(r) = J^^k'^F^R^ik'T) (40) 
10 



where R\ is the free particle solution, which can be determined from eq (I3ip by taking the 
limit x — » 0: 



i? (A=0) 



t l+l 



f(A) 



-iz/2 i\ 

T(aT 



r (1 + 1 + 2zio) 



r(2iio) 



2 F X (1 + Z, 1 + Z + 2iw; A; t) 



r (1 + I + 2iw) 



r(2wt>) (2w) 



i+i 



+ z 



(41) 



(42) 



It follows from the above that the density of states remains unchanged, as assumed in eq 
In the limit w —* oo, the correct centrifugal term is recovered and so it is expected that the 
standard radial Fourier transform, as given in eq (|18p. is also recovered. This result can be 
verified by noting the following identity: 



ii(*/2) 



-iz/2 i\ 

T(aT 



and considering the limiting behaviour of the T function and 2 F\ function, given in eq (J36J) 
and eq (|37|) . The standard transformation is indeed found in the limit w — > oo. 

The same derivation of the Sommerfeld factor as presented in Section [5] can be carried out 
up until eq (|17p , where the Sommerfeld factor is given as an integral of the momentum space 
wavefunction weighted by various powers of momentum. At this point, to evaluate the integral, 
it is necessary to consider derivatives of eq (|40p instead of the standard radial transformation. 
Using the series expansion in eq (|42f) for the free particle solution, the Sommerfeld factor for 
a two body system with the modified centrifugal potential, S l: is then: 

2 



s, 



(21 + 1)! T(2iw) (2w) 



' 1 d l Riiz) 



(I!) 2 T(l + l + 2iw) dz l 



z=0 



(43) 



Note for I = and/or w — > oo, this reduces to the usual Sommerfeld factor as expected when 
the centrifugal term and/or its substituted approximation vanishes. This can be seen using 
the relation T(l + h) = hT(h) for the I = case, and eq (I36D for the w — > oo limit. 



3.2 Analytic Sommerfeld factor 

The approximate analytic Sommerfeld factor is now presented for the case of Yukawa interac- 
tions for arbitrary partial waves. Substituting the approximate wavefunction, given in eq (|3ip . 
into the Sommerfeld factor which accounts for the modified centrifugal term, the following 
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result is found: 

r(g-)r(a+) 1 2 

r(l + I + 2iw) l\ 



In the Coulomb limit, this becomes: 



Using the relation T(l + h) = hTih) recursively, and that |T(1 + ib)\ = y/nb csch(7r6) for 
real b, the Coulomb case can equivalently be written as: 

Six* = S x fl (l + where S = (45) 

Recall that this result is exact for the Coulomb limit. Although the perturbative cross sec- 
tions of higher partial waves are suppressed by factors of (3 21 in the non-relavistic limit, this 
result indicates that in the zero velocity limit (x — * oo), the Sommerfeld effect introduces 
a multiplicative factor that leads to all partial wave cross sections having the same velocity 
dependence. 

In the large screening limit for the Yukawa potential (w — ► 0), the Sommerfeld effect 
saturates when the de Broglie wavelengths of the scattering particles are much greater than 
the range of the potential. However, in the attractive case, resonances can develop when the 
system is close to bound states leading to large enhancements, and these states can form at 
zero energy. For finite range potentials such as the Yukawa case, the number of bound states 
present is always finite. Below a critical value for the coupling, no bound states can form and 
so the magnitude of the Sommerfeld factor is bounded. In the limit w — > 0, the Sommerfeld 
factor is approximately: 

5, |r(i + z + v^) r(i + /-V^)| 2 /G0 4 (46) 

where to = y (m^/5). The Gamma function has poles at and the negative integers. There- 
fore, the above result suggests that resonances are found when wx = (1 + I + n) 2 with n 
being a non-negative integer. The critical values of y for the exact and approximated Yukawa 
potentials that allow one bound state of angular momentum / = 0, 1 are given in Table 1. The 
resonance behaviour at these points will be demonstrated in Section [H 




r [1+1+ 



IX 



(l\ 



H2 



(44) 
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Potential 


I = 


I = 1 


approx Yukawa 


5/ 


45/ m+ 


exact Yukawa 


1.680 


9.082 



Table 1: Critical values of y (= AM/m^) for one bound state with angular momentum I. In 
this paper, 5/m+ = 7r 2 /6 (~ 1.645) has been used. The exact results are taken from |16j . 

In the large screening limit, for small wx, the approximate analytic Sommerfeld factor 
behaves in a series expansion as: 

S t = l + 2wxVi(l + + 0(w,wx) 2 (47) 

where ip^ is the trigamma function, defined by: 

Mh) = ^lnT(h) (48) 

For physical angular momentum, ipi(l + 1) ~ 1/(1 + \) is a good approximation, exact in the 
large / limit, and 18% (3%) accurate for Z = 0(1). This result suggests that for a given I, the 
magnitude of the deviation of the Sommerfeld factor from unity is roughly equivalent for the 
attractive and repulsive cases in this limit. The same behaviour is also seen in eq (|45p for the 
small x Coulomb limit, independent of I. 



4 Comparison of numerical and analytic results 

In this section, the Sommerfeld factor determined from numerical simulations and using the 
approximate analytic result is presented and compared. The numerical results were found 
with the unapproximated Yukawa potential and centrifugal term. Fig 3 shows the Sommerfeld 
factor for the I = partial wave cross section. 

The Coulomb limit is reproduced in the limit y — > oo. For the attractive case, resonances 
are found with the critical points in agreement Table 1. In the repulsive case, there appears to 
be a qualitative symmetry between the large screening and Coulomb limits. This is expected 
by the similar forms that the approximated Sommerfeld factor takes in these limits. The 
transition between the two regions is roughly where the contours would intersect if the limiting 
results were extended. By equating the first order terms in a series expansion of eq (|45p 
and (|47|) . this is found at y ~ (21 + l)x. 



13 



10 



10 



1.3*- 



1.1 



10 -2 10 _l 1 10 10 2 10 3 

x (= A//3) 

(a) Attractive Potential 



lO" 1 - 



10 



10 



0.5 



lO" 1 1 



10 10 2 



10 3 



x (= A/y3) 
(b) Repulsive Potential 



Figure 3: Sommerfeld Factor for a Yukawa interaction of a s-wave state where the contours not 
labelled vary by a factor of 10. The solid lines represent S Q with the wavefunction determined 
from numerical simulations, and the red dashed lines are the approximate analytic result, S . 
These contours are almost coincident except around the resonances where there is a noticeable 
shift. For clarity the contours of S > 10 3 in (a) are not shown. 



The analytic factor is accurate to within 10% in the regions where S < 10 and S > 0.1, 
with the largest errors occurring when y/x ~ 1. This is not surprising as the choice of <5/ m * 
is only accurate in the limit y/x — > 0, and the exact Coulomb result is recovered as y/x — > oo. 
The resonances are slightly off-position, with worse correlation for the higher resonances. 

Figd]shows the Sommerfeld factor for the I = 1 partial wave cross section. The correlation 
is relatively poor in the small y/x region. For y/x <C 1 and y = 1 or 3, the error of the 
analytic factor in reproducing the numerical result is roughly 10% and 50% respectively. For 
the attractive case, the incorrect resonance positions for the approximated result lead to order 
of magnitude errors for y > 6. However, in the repulsive case, order of magnitude errors only 
are found for y > 30, with the analytic factor being always less than the numerical result. 
The approximate analytic factor is then of limited use when the Sommerfed effect is large, 
without finding a better approximation. 

If the shift 5 — * 1.2 5 is made, the contours at small y/x are nearly coincident with the nu- 
merical result. Actually, when the numerical results are determined for the Hulthen potential 
with an exact centrifugal term, the contours very closely match the Yukawa case without this 
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x (= A//?) x (= A//3) 

(a) Attractive Potential (b) Repulsive Potential 

Figure 4: Sommerfeld Factor for a Yukawa interaction of a p-wave state where the contours not 
labelled vary by a factor of 10. The solid lines represent S 1 with the wavefunction determined 
from numerical simulations, and the red dashed lines show the approximate analytic result, 
S 1 . The contours are almost coincident in the region y/x > 1 but have a mismatch in the 
region y/x < 1. For clarity the contours of S 1 > 10 4 in (a) are not shown. 

shift. The expected origin of the error in the analytic factor therefore is dominantly from the 
approximation applied to the centrifugal term. A criterion has not been found to better select 
5 once this approximation is applied, but the following guess at a substitution was determined 
to significantly improve the results in the I = 1 case: 

w — ► w [1 — c x tanh (c 2 /w)] 

For the choice (c l5 c 2 ) = (1/6, 1/9), the results then agree to within 5% in the regions where 
Si < 10 and Si > 0.1. The largest fractional error for the approximated analytic result in 
these regions still occur at y/x ~ 1. This transformation has effectively shifted the y variable 
in the small y/x region and leaves it unchanged in the large y/x region. The position of 
the resonances is still slightly offset though leading to large order of magnitude differences 
between the results where these features are present. 

For the attractive Yukawa potential, the resonances of different partial waves do not co- 
incide. It is then possible for a system to be dominated by a higher partial wave if close 
to a large I resonance. For the repulsive case in the large screening region, the Sommerfeld 
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suppression is reduced for higher partial waves since the centrifugal barrier becomes more 
effective at keeping the wavefunction away from the interaction core. It is found using eq (I46p 
that the reduced Sommerfeld suppression for a successive partial wave compensates for the 
perturbative f3 21 suppresion, (3 2 {S l /S l _ 1 ) > 1, when: 



where some 0(1) factors have been neglected, and the second inequality follows since the 
result only applies in the large screening region, y < (2Z + 1) x. To determine for a particular 
case how many partial waves are significant and which channel dominates, it is necessary to 
also consider the relative magnitudes of the velocity independent factors in the partial wave 
expansion. For an attractive Yukawa interaction in the small y region, the larger angular 
momentum processes are enhanced less by the same argument, so higher partial wave terms 
which are negligible in the perturbative expansion stay negligible. 

As is suggested by the numerical simulations, and also by the approximate analytic results, 
the Sommerfeld factor tends to infinity as /3 — ► 0, when approaching resonances of the system. 
In physical situations, this would not occur as the zero energy bound states generally have 
finite lifetimes. A more careful analysis is then required in the close neighbourhood of these 
resonances. 

5 Conclusions 

The Sommerfeld effect can bring large order of magnitude corrections to cross sections. In this 
paper, the Sommerfeld factor relating the non-perturbative and perturbative matrix elements 
has been derived for arbitrary partial wave channels. It is demonstrated that in the non- 
relativistic limit, the incoming states are transformed by scattering interactions to effectively 
form a composite state whose wavefunction is a solution of the relevant Schrodinger equation. 
This wavefunction weights the perturbative matrix element in an integration over momentum 
space to give the non-perturbative matrix element. 

The Sommerfeld factor for Yukawa interactions has been determined in an approximate 
analytic form for arbitrary partial waves, and also evaluated by numerical simulations for 
the I = 0, 1 cases. The s-wave result is found to be accurate to within 10% when the non- 
perturbative cross section is up to an order of magnitude different to the perturbative result. 



l-[3 



< 




(49) 
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For the p-wave factor, a correction introduced to compensate for the centrifugal approximation 
allowed the analytic result to be accurate to within 5% for the same case. This is then a 
promising approach for application of the Sommerfeld effect when computational resources 
are prohibitively restrictive. However, the analytic result is poor at reproducing the correct 
resonance structure for attractive potentials. If the resonance region is probed, it remains 
necessary to use numerical simulations in order to obtain accurate results. The results are 
exact though in the Coulomb limit. It was found for certain areas in the parameter space that 
higher partial waves can dominate cross sections, when a perturbative analysis suggests that 
they are negligible. It can therefore be critically important to include the Sommerfeld effect 
for the I / case. 

One consequence of the Sommerfeld effect is the possibility that the annihilation mech- 
anism which controls dark matter freeze-out could be different to the dominant mechanism 
of present day annihilation in our galactic neighbourhood. This is of relevance to identifying 
potential signals for indirect dark matter detection. For this case, "boost" factors are com- 
monly required in order to find a detectable signal. The Sommerfeld factor of an attractive 
interaction could provide such an enhancement. This can also be present in addition to other 
mechanisms such as internal Brehmstrahlung events or dark matter density perturbations. 

During preparation of this paper, a similar derivation of the general Sommerfeld factor 
and evaluation for the Coulomb case has been presented in [10] , in agreement with the results 
in this paper. In [11] . the s- and p-wave Sommerfeld factors were also evaluated numerically 
given an attractive Yukawa interaction for various slices in the parameter space. Those results 
are also in agreement with that presented in this paper. 
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